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We construct the scalar and graviton propagator in quasi de Sitter space up to 
qh| first order in the slow roll parameter e = —H/H^. After a rescaling, the propagators 

are similar to those in de Sitter space with an e correction to the effective mass. The 
fDJ), limit e — > corresponds to the E{3) vacuum that breaks de Sitter symmetry, but does 

not break spatial isotropy and homogeneity. The new propagators allow for a 
self-consistent, dynamical study of quantum back-reaction effects during inflation. 

> 

a^ ■ 1 Introduction 

cn 

■ While inflation -the accelerated expansion of the early universe- was originally 

! introduced to solve the horizon, flatness and cosmic relic problem [T][2], it was soon 

I realized that inflationary space-times have the interesting property of particle 

creation [3] [5]. This can be seen from a simple application of the energy-time 
uncertainty relation [6] [7] that states that a virtual particle-anti particle pair with 
^ I energy E and momenta ±fc can exist a time 6t given by 

^ -t+5t 

dt'E{t\k)<l. (1) 



While in flat space 5t is very small, it can grow large in curved space-time. In 
particular in the context of cosmological inflation, one finds that for massless modes 
with wavelengths larger then the bubble radius, 5t can become infinite. This means 
that during inflation massless, infrared particles are created out of the vacuum. Since 
photons and massless fermions couple conformally, their production rate is suppressed 
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by the scale factor. The only non-conformally coupled, massless particles are massless 
minimally coupled scalar s and gravitons. 

Based on the discussion above one must conclude that quantum effects from 
massless minimally coupled scalars and gravitons during inflation can be very 
substantial. Indeed, these effects are exactly the source for scalar and tensor 
cosmological perturbations [3]. The imprint scalar perturbations leave on the cosmic 
microwave background has been verified by WMAP [8]. Tensor perturbations might 
be detected by the Planck satellite. 

In general doing quantum field theory on a curved background is very involving, 
but, due to its high amount of symmetry, much more can be done in the context of de 
Sitter space (section [2]) . Since in de Sitter space the scale factor grows exponentially 
with time, it is an excellent paradigm for an infiationary space. The L'-dimensional 
propagators in de Sitter space are known|9j tlOjpJJp^ and therefore it is possible to 
use dimensional regularization to calculate many kinds of interesting and possibly 
observable quantum effects during inflation. For example a massless minimally 
coupled scalar with a quartic self-interaction has been studied [10] [13] |14j . The 
resulting model shows a violation of the weak energy condition on cosmological 
scales. A non-minimally coupled, massive scalar with quartic self-interaction has also 
been studied [T3]. In this model the radiative corrections to slow roll inflation where 
calculated and found to be unobservable. 

Also scalar electrodynamics has been studied extensivelv[6] [T6] [T7] [18] [T9] [20] . The 
vacuum polarization has been studied and it has been shown that in such a model 
superhubble photons acquire a mass. The contribution to the zero-point energy of 
these photons sources cosmological magnetic flelds. Most recently the theory has 
been formulated stochastically [21j. 

Another interesting line of research is the study of quantum back-reaction on the 
vacuum[22]p3]pl|[25]p6][27][28]. In a pioneering study [22] [23], it was shown that 
the two loop contribution to the graviton one-point function slows inflation. The 
resulting model, where inflation is driven by a large cosmological constant and ends 
naturally due to quantum effects is an alternative to scalar fleld driven inflation. 
Moreover it might prove to be important for a better understanding of the 
cosmological constant problem [29j. However this calculation was done in 3 -|- 1 
dimensions, using a (time dependent) ultraviolet momentum cutoff. Since it is not 
clear whether the observed effect is a relic of the regularization procedure, it is 
important that this calculation is redone using dimensional regularization [30] . 

There are however several drawbacks to the use of de Sitter space. First of all, 
while de Sitter space provides an excellent framework to study properties of inflation, 
it is never truly realized in nature. Indeed, since in de Sitter space the Hubble 
parameter is per definition globally constant, it follows that if the universe was once 
de Sitter it is always de Sitter. From the fact that the universe is not de Sitter today, 
it follows that it was never de Sitter. The deviation of the infiationary universe from 
de Sitter space is even measurable through the deviation of the spectral index 
from scale invar iance [2] [8]. When these deviations are small enough, one could work 
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in a space that is 'locally' de Sitter, however in that case there is no control on the 
error one is making. 

The second problem with de Sitter space is that it is non-dynamical. While 
space-time is certainly curved, all dynamics are governed by the Hubble parameter 
and the Hubble parameter is constant. Strictly speaking it is therefore inconsistent to 
study how space-time changes due to back-reaction effects if the calculation of these 
effects explicitly assumes that the background is constant. To self-consistently 
consider back-reaction, one must allow the Hubble parameter to change with time. 
Although in [22j ^23j a distinction is made between the constant Hubble parameter H 
and the - dynamical- physical observable -ffeffj it is still not fully consistent. The 
reason is that the propagators are calculated using the constant H, where they 
should be calculated with the dynamical, time dependent -ffeff- 

The final issue with de Sitter space is the fact that it is not possible to construct a 
propagator for a massless minimally coupled scalar or for the graviton that respects 
the de Sitter symmetry |31j. Usually this problem is circumvented by letting these 
propagators break the de Sitter symmetry, but not spatial isotropy and homogeneity. 
The argument is that it is these symmetries one expects to be present in the early 
universe, not the full de Sitter group (because of the problem described above). 
However, if we are forced to abandon de Sitter symmetry anyway, why not do the 
complete calculation in a space with only spatial isotropy and homogeneity, but 
otherwise similar to de Sitter space. 

Constructing the propagator in this space, known as quasi de Sitter space, is 
exactly the aim of this paper. Quasi de Sitter space is a generalization of de Sitter 
space, where H depends only mildly on time (see section [2T]) . While this time 
dependence breaks the full de Sitter symmetry, it preserves spatial homogeneity and 
isotropy. We therefore find that there is no problem to describe massless minimally 
coupled scalars or gravitons in quasi de Sitter space. Moreover, since quasi de Sitter 
space is both dynamical and a more realistic approximation in many inflationary 
models then de Sitter space, the other two problems mentioned above are also 
resolved. 

In section [2] we overview some properties of de Sitter and quasi de Sitter 
space-times. In sections [3] and H] we calculate the scalar and graviton propagators 
respectively. In section H] we also calculate the ghost propagator associated with our 
gauge fixing. We conclude in section [5l 



2 de Sitter and quasi de Sitter space 

4-dimensional de Sitter space is the hypersurface given by 

-Xi + Xf + Xl + Xi + Xl = ^ ( 2 ) 

embedded in 5-dimensional Minkowski space-time, where H is a constant parameter 
known as the Hubble parameter, de Sitter space is a solution to the Einstein 
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equations with a positive cosmological constant, in which case 

= f . (3) 

The isometry group of de Sitter space, 50(1,4), is manifest in this embedding. We 
win use flat coordinates, which cover only half of the de Sitter manifold, given by 
(i = 1,2,3) 

ri L 

1 H 

X4 = — cosh{Ht) - —Xix'e"\ 
H 2 

—00 < t,Xi < 00. 

In these coordinates the metric reads 

5^,, = diag(-l,a^,a^,a^) a = (5) 

such that a/ a = H. The coordinates ^ can be written in conformal form by 
changing coordinates to conformal time r] defined as adr] = dt. The metric becomes: 

g,iu = a^Vfiu, ^/ii^ = diag(-l, 1,1,1), (f? < 0) . (6) 

We define the de Sitter invariant distance functions 

Z{X;X)=H^7]abX^X'' = 1-^Y{X;X). (7) 
In conformal coordinates these functions read 



1 

2' 

Axl+{x;x) (8) 



z++{x;x) = i - -y++{x;x 



y++[x;x = 

rjrj 

/S.x\^{x; x) = —{\r] — fi\ — ie)"^ + \ \x — x\\'^ , 

where Y{X;X) = y(x;x), Z{X;X) = z{x;x), are the functions given in ([7]). a = a{r]), 
d = a{fi) and e > refers to the Feynman (time ordered) pole prescription. In order 
to apply the Schwinger-Keldysh (see e.g. the appendix of [32]) formalism we also 
define 

Ax^ {x;x) = —{rj — f] + ie)"^ + \\x — x\^ 

Ax^^(x; x) = —{rj — fj — ief + \\x — x\^ (9) 
Ax^_(x; x) = —{\ri — vi\ + ie)^ + | |x — x| |^, 
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from which , y |- and y fohow immediately. For the rest of this paper we use 

y = = A4+ (10) 

unless stated otherwise. When e = 0, the function y is related to the geodesic length 
I = i{x; x) between points x and x as, 

y{x;x) = Asin^ ^H£{x;x)^ . (11) 

The function y has the following properties: if y = 0, 5; lies on the lightcone of x. For 
y < 0, x lies in the future or past of x and for y > 0, x and x are space-like separated. 
If y = 4, X lies on the light cone of the antipodal poin10 of x and if y > 4, x is in the 
future or past of the antipodal point. 



2.1 Quasi de Sitter space 

Quasi de Sitter space is defined by the requirement that H depends only mildly on 
time: 

H = Ho-Hlet + 0{e^), e ^ _^ « i, e = ©(e^). (12) 

While such a Hubble parameter does not correspond to an exact solution of 
Einstein's equations, it is a good approximation in most inflationary models. For 
example in scalar field models, where inflation is achieved by a minimally coupled 
scalar field 4>, slowly rolling in a potential V, e is just the slow roll parameter 

1 ( d{y{ct>))/d<t> ^^ 

and (|12|) is correct up to order [2] . The metric of quasi de Sitter space can be 
written in conformal coordinates as 

9^ll^ = a^r]f,^ (14) 

where the scale factor is given by 

a = (l + (1 + Ha))e) + 0{^) (15) 

Since we know from ^ that t = ln{a)/HQ + 0(e), this is equal to 

a = -^{l + e) + 0{e') (16) 

If we take flat spacial slicings, the symmetry group of quasi de Sitter space is E{3), 
which corresponds to spatial homogeneity and isotropy (invariance under translations 
and rotations). 



^The antipodal point is obtained by the transformation of the coordinates (U) Xa —Xa, with a 
(0...4) 
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3 The scalar propagator 



The Feynmann propagator 



(17) 



of a minimally couplecH scalar field with 
equation: 



mass m satisfies the following Klein-Gordon 



^/—g{\I\ — m'^)iA{x; x) = i6^{x — x) 



(18) 



where □ denotes the d'Alembertain 



a 



(19) 



and 5^ is the /^-dimensional Dirac delta. In de Sitter space, one can solve for A(x;x) 
purely in terms of the de Sitter invariant distance function y defined in dH]) [9] [31] . 

Replacing with y gives the anti-time ordered propagator (which has the 

opposite sign in front of the delta function) , while replacing it with y_| and y 

gives the two Wightman propagators. The result, known as the Chernikov-Tagirov 
propagator 0, however becomes singular in the limit when the mass goes to zero. It 
was shown in [31][33j that it is not possible to construct a physically meaningful de 
Sitter invariant propagator for the massless, minimally coupled scalar field. The 
reason is that in such a case the solution to (llSp has an additional singularity when x 
is on the lightcone of the antipodal point of x. Such a solution corresponds to the so 
called a-vacuum. Although the issue is disputed, these vacua are usually considered 
unphysical [3l!|35][36] [37] [38]. 

In quasi de Sitter we find one can do something similar. It turns out that after a 
simple rescaling, we can construct a propagator purely in terms of the function y. Due 
to e corrections, this propagator now has a consistent massless limit. The y function 
we use is the one defined in ([8]). One could think however of different definitions for 
y. For example the most general definition appears to be through the geodesic length 
(jlip . However the geodesic length in quasi de Sitter turns out to be a rather complex 
object and not very helpful in solving (|18|) . An alternative definition of y is 



While this definition is equivalent to dHj) in de Sitter space, it differs by a factor 
(1 -|- 2e) in quasi de Sitter. We use the definition ([8]), since in that case the 
singularities of the general solution to (|18p are at y = and y = 4 and thus resembles 
the de Sitter solution most. 

Our rescaling of the propagator is 



^Our results straightforwardly generalize to a non-minimally coupled scalar by replacing m? — > 
with R the Ricci scalar and ^ a a coupling parameter. 



{x;x) 



(20) 



A{x;S:)^{adrE{y) 



(21) 
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where it; is a constant parameter. We shall determine w by the requirement that H is 
a function of y only (strictly speaking H is not a function, but a distribution). We 
rewrite the left hand side of (jlSp in terms of y-derivatives: 



'dy' 



+ 



where 

and we used 



diy + dfy) -[D-2 + 2ew) (-) {doy) 



I? = {D- l)ewHl + + 0{e^) 



d 2 2 

— -a ^ 
dy 



(22) 



H + 2H\ 



(23) 



(24) 



An important point of this calculation is that, even if the mass of the field is zero, 
the effective mass jj. as it appears in (j22p is nonzero. Therefore we find the propagator 
of a massless scalar field in quasi de Sitter space does not suffer from the presence of 
a-vacua. Indeed this is expected, since in order to construct the propagator for such 
a field in de Sitter space, one exactly breaks the de Sitter symmetry to the symmetry 
of quasi de Sitter space. The several terms in (f22|) evaluate to 



-2 



- {doyf + {d,yf 
dly + d'^y) -{D-2 + 2ew) (-) {doy) 



-_H\l-2e)y{4-y) + Oie^) 
-.-H\l-2€)(D{y-2) 



(25) 



+ e{D-2 + 2w)y 



where we made use of (I15p and (I16p . The (5- function is only sourced by the singular 
term in H for y — > 0. Putting everything together we find for the nonsingular terms 
that 



y(4-y)( 



d_, 
dy' 



{D{y -2) + {D-2 + 2w)ye] 



d 
dy 



[l + 2e 



if2 



^H(y) = 0. (26) 



In order for H to be a function of y only, the left hand side of (j26[) needs to be a 
function of y only. The term fj? /H'^ can in general never be written as a function of 
y. However if we assume that our field is light, such that 



m 



0{e), 



(27) 



' / H"^ is, up to order corrections, proportional to fi^ /Hq and thus a constant. The 
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resulting equation is a hyper geometric equation with a general solution 

, , ^rD-l D-2 + 2W D-1 D-2 + 2W Dy. 
ir.{y) = A2Fi[-^^ + ^ e + z^D,^— + ^ ^"^^'^'4^ 

^ ^ rD-l D-2 + 2W 
+ S2Fi[^— + e + i/D, (28) 

D-ID-2 + 2W D , , y, 

+ ^ e-un;- + {D-2 + 2w)e;l-^], 

where 

2 1 Z)-2 + 2?i; \2 o A /"^ 



The constants A and B can be fixed by looking at the singular terms. First of all we 
require that the propagator has no singularity at y = 4. Such a singularity would 
correspond to an additional delta-source at the antipodal point (a-vacuum) and is 
considered unphysical. This requirement fixes ^ = 0. We fix B by requiring that the 
singularity at y = sources the (5-function correctly. The 5 function is sourced by the 
most singular term, which is given by 

^ _ r.fyy-''^^ r[D/2 + {D-2 + 2w)e]T[D/2 - 1] 

The delta function is sourced by letting ^/—gO act on this term. We find 

z6^{x -x) = (aa)^"'a^-2a2zH(y)sing 

_ (a)^-"r[Z)/2 + (Z)-2 + 2^;)6]r[Z)/2-l] .d-2^2( 1 ^ 

(31) 

where Ax^ is defined in ([H]). Depending on the propagator we calculate, we use one of 
the following relations 

If one uses the appropriate form of ()18|) they all give the following result 

T[D/2 + {D -2 + 2w)e] (47r)^/2(l + e)i5-2a2e«, • ^ > 



We can now fix w in order that B is independent of a and rj. To do this we write 

rD-2 



^ = H^-\l -e{D-2 + 2w) ln(a)) + 0(6^) (34) 



and therefore if we choose 



w 



D 

1 

2 



(35) 



H is truly only a function of y up to order e'^ corrections. Remarkably, for this choice 
of w, equations (f26j) and (f28]l reduce exactly to their analogues in de Sitter space, 
with only the mass term modified. The final propagator is 



iA{x; x) 



D-2 prP-1 I ,,„irr£'-i 



(1 + 

X 2F1 



\D-2 



D 



(4^)^/2 

1 D 
— + — ; 



r[f] 

D 



(36) 



with 



-I) 



D 



1\2 m2 

- 772 + o(^-l)(^-2), 



m 



0(e) 



(37) 



The A , and A |_ propagators can be obtained by replacing in ([36]) by 

y , y_| and y |_, respectively. The reason is that the constant B, given in (f33l) . is 

the same for all propagators. 



3.1 The de Sitter limit 

In the limit e — > 0, the leading order term for the massive scalar propagator is 



lim iA(x: x) — , , „ ,^ 
e^o ^ ' ^ (47r)^/2 

X 2F1 



D-l D-1 D y 



(38) 



+ Oie) 



with 



lim Vri 



D-l\2 m2 



Hi 



(39) 



This is ~ as expected - the Chernikov-Tagirov propagator [9]. For a massless scalar 
the story is slightly more complicated, since the term 



1 

r[^^ vd] 



(40) 
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becomes singular. We split off this singular term by rewriting the hypergeometric 
function in (f36j) as a series: 

(aa)-^(^/2-i) H^-^ T[D/2]T[l-D/2] 



A(x, x) 



l + e)^-2 (4^)^/2 r[l/2 + z/fl]r[l/2-z.i5 

n=l 

oo 



n= 

r[ 



T[D/2 + n]T[n + l] \A. 



n=-l 

r-D-1 , ,. i-nrD-l 



r[3-L>/2 + n]r[n + 2] V4 



(41) 



T[D/2] 

Now we can easily expand around e = and obtain 



m=0 



D-2 




(47r)^/2 

r[D - 



-2r[D - 1] 1 



+ 



1] 



+ E 



T[D/2] 
T[n + D 



- 2)r[Z)/2] e 

In(ad) + 2-JE + '4'{D - 1) + vr cot 



(42) 



^ nr[I?/2 + n] V4. 

r[n + D/2+l] 



(y_ 

^ {n- D/2 + 2)T[n + 2] \A 



ri-D/2+2 



+ 0{ 



where 7^; is the Euler-Mascheroni constant and il)[z\ = ^ ln(r[z]) denotes the 
digamma function. Apart from the finite constants 2 — 7^; + 'iIj{D — 1) and the 1/e 
term (which is also a constant in our approximation), ()42|) exactly corresponds to the 
massless scalar propagator in [lO]. Since a constant corresponds to a solution of the 
homogeneous equation (scalar condensate), it can always be added. This can be seen 
from the fact that in the massless case the hypergeometric function multiplying in 
(f28l) is to leading order just 1 and thus we are not forced to put A = Q. The 
propagator in [lU] is calculated by requiring that the inevitable breaking of the de 
Sitter symmetry (to avoid a-vacua) by a massless scalar, is such that it does not 
break spacial homogeneity and isotropy (this is known as the -^(3) vacuum). Since 
those symmetries are the symmetries of quasi de Sitter space, it is expected that the 
de Sitter limit of the quasi de Sitter propagator would correspond up to a constant to 



4 The graviton propagator 

Various attempts have been made in the past to calculate the graviton propagator in 
de Sitter space. In many works [39] [lO] [41j |42j this is done by adding a de Sitter 
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invariant gauge fixing term. However, it turns out that these Green's functions give a 
divergent response to a point mass and moreover they do not solve the classical 
invariant equation of motion [13]. Based on this Tsamis and Woodard conclude that 
there must be some deep, still not completely understood physical principle that 
forbids the use of de Sitter invariant gauge fixing terms [llj^44j [45j. This problem 
might be related to the issues concerning the linearization instability in de Sitter 
space [TT] [44 j |47j [48 j [49 j . The point is that the linearization instability implies a 
constraint on the (classical) field equations. This constraint has a nonzero response to 
sources on the whole de Sitter manifold. However, if one then adds a de Sitter 
invariant gauge fixing term and calculates Green's functions, one finds that these 
Green's functions only respond to sources in the past light cone. Since this does not 
cover the whole manifold, one has changed classical physicsH 

This problem is circumvented by working only on half of the de Sitter manifold 
and by using a gauge fixing term that breaks the de Sitter symmetry [H]. If one 
calculates the graviton propagator in this way, one finds that the solution can be 
written in terms of three scalar propagators, multiplying some tensor structure 
^30j^44j. Since one of these propagators corresponds to a massless minimally coupled 
scalar, the graviton propagator inevitably breaks the de Sitter symmetry. 

In this section we calculate the graviton propagator in quasi de Sitter space. 
Although slightly more complicated, the result is similar as in de Sitter space. Two 
differences are that our gauge fixing term, the analogue of the one used in [30j [44j . 
does not break the quasi de Sitter symmetry. Moreover, as we have shown above, a 
massless minimally coupled scalar poses no problems in quasi de Sitter space. 

4.1 The graviton action 

We start our calculation with the Einstein-Hilbert action with a cosmological 
constant A, coupled to a scalar field </>. 



and the hats indicate 'background+perturbation'. We have added the matter action 
Sm for consistency, since without it, the quasi de Sitter metric is not a solution to 
Einstein's equations. Notice that the addition of matter induces a mixing between 
the graviton and the scalar field. Taking this mixing into account is crucial for the 

■^The problem seems to arise as soon as one imposes an average gauge condition, as when one adds a 
gauge fixing term. In Proca theory the de Sitter invariant propagator has been calculated using an exact 
gauge condition [12]. This propagator does give the correct response. 




(43) 



where 




(44) 
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calculation of scalar cosmological perturbations ^50j. The four dimensional 
propagators, with mixing, are calculated in [51] [52] in power law expansion. A study 
of the mixing between the gravitational and the matter section in D dimensions will 
be done elsewhere ^53j. Quasi de Sitter is a good approximation if the scalar field (p is 
slowly rolling in the potential V [2j. Using the geometric quantities from appendix [Aj 
the Priedmann and fluid equation corresponding to are found to be 



(45) 



where a dot denotes a time derivative. By choosing the appropriate potential, ()45p 
can describe any time evolution of the scale factor. The action (|43p therefore 
describes - as far as the geometry is concerned - a completely general FLRW 
universe. To calculate the graviton action we decompose our metric in a background 
g plus a perturbation h [M] : 

If we require that g^'^gfip = 5'^ and we raise and lower indices on the perturbation 
with the background metric: h^a = we find that 

6g^"' = -h'"' + /I'^^/i"^ + 0{h^). (47) 

We expand the action in powers of the perturbation around the background as 

S = 5o + 5i + ^52 + ..., (48) 

where it should be noted that, since 5g^^'^ is quadratic in /i. Si contains a part 
quadratic in h. For the first two terms we get 



So = ^jV^[R-{D- 2)a) + Sm 



(49) 



where 



2 6Sm 
^^9 



fJ,U 



{EE)^iy = R^u — -g^uR H ^ — ^9fiu — i^^^t^y 



{d^<t^){dA) - 9J\{dact>){dp<l^)9'''' + y (</.)) 

^ (50) 



12 



If we drop terms cubic in h, S2 is given by 



So 



(51) 



where we have defined the trace reversed graviton: 
Using the following identities [55] 

and removing boundary terms, we find 

1 



Q/3 



(52) 

(53) 
(54) 



52 ^1 

2 ~K 



-9 



(55) 



where the d'Alembertian is given in (I19p . Adding the part quadratic in h from Si the 
total quadratic part of the action becomes: 



1 



h'''' ( Bg^^gp, + 2R^pp, ) h''" + l(V^/i^")(V^/i, 



Ir, 



+ |(/i'^^(5V) - Tp^h'^^K^^) + -hP''hp>^{EE)p, 



(56) 



In the following we drop the {EE)py terms from (I56p . since they are zero on-shell. 
The contribution to the graviton action from the matter action can be found from 
(IHnll and 



5T, 



to give 



K 



KA\idc.(t>){dp(t>)9''^ + V{ct>)) + ]-gp,{d^cl>)[dpcl>)h'^'' (57) 



Prom (1451) we can derive that 



K{(j)'f = -2{D - 2)Ha' 



(58) 
(59) 
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where a prime denotes a derivative with respect to conformal time. Using this, the 
graviton action becomes 



5p) = - 



+ 



D-2 . 



H(hh°o 



(60) 



In any FLRW metric we can choose coordinates where the metric is conformally flat, 



(61) 



A lengthy calculation with several partial integrations and using the quantities given 
in Appendix Ogives apart from boundary terms the following result ^44j : 



S{2) = „D+4 



iv'^-(5-5„-(D-2)(^)9o)V'M. 



L> - 2 / / a' \ 2 



D-2 



2 -T --o+HU^tl; 



(62) 



where 



We follow [44J and add the following gauge fixing term: 



a^+4(V^^'^°)(V,V^^ 



(63) 



(64) 



While this gauge fixing term breaks the de Sitter symmetry, it does respect the quasi 
de Sitter symmetry. Removing another total derivative and, using ([MI)) we obtain the 
general form of the gauge fixed graviton action in FLRW space-time 



S. 



(2) 
GF 



1 



D-2 



{H^ + H)5l5P5"A^,^, (65) 



where 



□ 



s 



a Oaa rj '^Op 



(66) 



is the d'Alembertian as it acts on a scalar field. 
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4.1.1 The Ghost lagrangian 

The ghost Lagrangian corresponding to our gauge fixing (I64p is given by 

/:ghost = -a^'V^'SF^ (67) 
where V is the anti-ghost field and 

(68) 



To calculate we consider an infinitesimal coordinate transformation (indicated by 



y^^' = yl" + kV . (69) 
Under this transformation the metric g (see (I46p ) transforms as 

9U(y) = aAv) - ^^{gc.u{y)d^V^ + g^M^uV^ + y'dpg^^uiv)) + 0(y2) (70) 
and thus we find for the variation of the rescaled graviton field 



(71) 



Since we let be our anti-commuting ghost field, the second line will be cubic in 
the fields, when combined with (|67p . while for the purpose of this paper we are only 
interested in quadratic contributions. So dropping the second line, we obtain 



6F^ = - g^^d^d.V' + {D - 2)9^^^)30^ - {D - 2)( ^ j^o.^" 



a 



(72) 



Now using (j24p and (|66p we find from (j67p 

/^ghost = a^r/,^y^(5^n5 - (I) - 2) (7^2 _ ij)<505-)y^ (73) 



4.2 Calculating the propagator 

It is a well known fact that the propagator should be invariant under the interchange 
of both legs 

iA{x;x) = iA{x;x). (74) 

Looking at (I65p we therefore see that, since □s' does not commute with it acts 
on the other leg as 

□5(r?"^/>a/3) = Osia'm (75) 
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Thus from (j65p we find that the gauge fixed kinetic operator is 



1 



pa 



(76) 



where, because of ()75p indices are effectively raised and lowered with r/"^. The 
Feynmann propagator is defined by 

P^/'^[zp,A°^](x;x) = ^<5(°5^)5^(x - x). 

Also indices of the propagator should be raised and lowered with the Minkowski 
metric r]^u. We solve this equation by making the following ansatz, 

zUA"'3](x,x) =a(x,x)5("5^) + h{x,x)5l/;^^5l^ + c(x, x)r?p,7?"^ 

+ d{x, x) (syX^ + 5o"<5o^p<x) + e{x, 3:)8yj^5l 



(77) 



(78) 



This is the most general form that [po-A"'^] can take, consistent with the 
symmetries [44J. We insert this ansatz in (|77p and, combining terms that multiply the 
same tensor structure, we obtain 6 differential equations (in the following a, 6, c, d and 
e are functions of x and x) 



"fl "v 



ins(a) = ia~^6^{x-x') 
^Ds{b) -{D- 2){H^ + H)ia + ^b) = 







\nsid) + {D- 2){H^ +H){c-d) = 



Us{'-^d+\h + \e 







0. 



' 4 

ins(e) - {D-2)H^^b-d + e 
To solve this system first of all note that ()8ip and ()83p just mean that 
d = -{D -2)c- a 

e = -{D-2)d-b = {D- 2fx + {D - 2)a - b 

So the only independent functions are a, b and c. We redefine them as 

a = 2Ao 

6 = 4(Ai-Ao) 

2 .(A2 - (Z) - 2)Ao), 



(79) 
(80) 
(81) 
(82) 
(83) 
(84) 

(85) 



(86) 



{D-3)iD-2) 

where Ao,i,2 are functions of x and x. We find that the system (fTOl - ISljl is solved by 
the three functions A„ obevingpT] 



(^-n(D-n-l)(l-e))A, 



z5^{x-x') 



(87) 
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where e is defined in (fT2]) . Equation ([87|) is just the scalar propagator equation (fTSj) . 
So, just hke in de Sitter space [Hj, we find that we can write the graviton propagator 
in any FLRW spacetime in terms of three scalar propagators, with different couplings 
to the background curvature. Inserting ()86p in the ansatz ()78p . we obtain the 
graviton propagator in any FLRW background 

2 



Ipcr 



where 



+ 



(r?,. + {D- 2)6X){v»P + {D- 2)6^,6, , 



{D-2){D-3) 



??M.' = r?M^' + <^X- (89) 
How to solve ()87p in a general FLRW space is not knowr0, but from section [3] we 
know how to solve these type of equations in quasi de Sitter space. In that case the 
functions A„ are given by ()36p with 

2 _{D-2n-lf 



'{D,n) 



+ e\^-{D - 1){D - 2) - n{D -n-l] 



(90) 



Notice that in calculating (f90ll we did not use (i37|) . since the requirement = 0(e) 
is not met. However, there is in this case no problem in using (|29p to calculate I'd, 
since for the propagators (fSHI) fJ.'^/H'^ is a constant. The graviton propagator (f88l) is 
the main result of this paper. 



4.2.1 Comparison to mode function analysis 

It is an interesting exercise to see how our results compare to the standard treatment 
of gravitons in terms of mode functions. For this we only consider the 'physical' 
graviton ipij, defined by requiring that 

V'om = 0, V = 0. (91) 
We see that in this case the action ()65p ieads to the following equation of motion 



d^-iD-2)-doU 



We rescale ^ to the conformally coupled field tp'^ defined as 
The equation of motion for tp'^ is in Fourier space 







^2 ^ {D-2)iD-i) ^a'Y ^ D-2 a" 



\ a 



a 



(92) 
(93) 

(94) 



^It is however not hard to solve them if the equation of state parameter w i 
in those cases, even though e does not have to be small, it is constant. 



is a constant, because 
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Using (jlSp we find that 

{^y=^, ^=^, (95) 



l + 2e g" _ 2 + 3e 
77^ o rj 

where t] is conformal time. The equation of motion reduces [50] [56j [57j to 



, _ (i?(^Lz^ , '°-^f-^> .).i . 0. (90) 



This equation can be solved and the result is 

i^'i^ik) = ci^Hi^\kri) + C2^HP{kr]), (97) 

where H^^''^^ are Hankel functions of the first and second kind respectively and ci^2 
are constants. The index 1/ is given by 

2 (D-lf (D-l)(D-2) , , 

y'= ^ 4 +^ ^- (98) 

Comparing this result with (f90|) . we see that is exactly equal to 2^(_d,o)- The reason 
that we recover the n = index, can be seen from ([88]) by the fact that only Aq acts 
on the physical graviton ^ij. This indicates that our propagator correctly reproduces 
the standard graviton two-point correlations, obtained from the mode functions. 

From this calculation, we can also find a relation between v and n^,, the spectral 
index of graviton perturbations 

d In P„ , , 

The graviton power-spectrum Vg is, apart from some normalization 

oc A;^-Vli(^)P- (100) 

We can calculate the power-spectrum in the infrared by expanding ()97p for small k. 
The result is 

i^lj{k) OC 7]^/\7]k)-'' + 0{k-^+^) (101) 

and therefore we find that 

ng = {D -I) -2u = -{D -2)e + 0{e^), (102) 
which in D=4 indeed corresponds to the well known result that Ug = — 2e [2] [50j. 

4.2.2 Ghost propagator 

By similar arguments as in the previous section, we raise and lower indices on the 
ghost field with rj^^'^ . It follows from ()73p that the kinetic operator is given by 

= a^is^^Us -{D- 2){H' - H)6^^6^). (103) 
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The propagator is defined by 



V^^i[^^P]{x- x) = i5P,5^{x - i) (104) 

and we find that 

iU^P]{x- i) = z5SAo(x; x) + z5°5^Ai(x; x), (105) 
where the A„ propagators satisfy 

[j^-^iD-n- 1)(1 + 6)) A„ = ^^^^ (106) 

Just as (|87|) . the functions A can be computed as in section [3l The solution is given 
by dMI) with 

'>(D,n) = ~ ^4 ~ + " ^^^^ " " ^^^^ (107) 

5 Conclusion 

While de Sitter space is an excellent framework to do calculations of quantum effects 
during inflation, it has two drawbacks. The first is that it is never realized in nature 
and the second is that it is not dynamical in the sense that the Hubble parameter is 
constant. To solve these problems one needs to do calculations in quasi de Sitter 
space, where first order correction in the slow roll parameter e = —H/H^ to the 
background geometry are taken into account. Not only is quasi de Sitter space a 
more realistic approximation of the inflationary universe, it is also a dynamical space. 
Indeed, since H is nonzero, the Hubble parameter is allowed to change. Especially if 
one wants to consider back-reaction effects or calculate quantum corrections to tensor 
cosmological perturbations, taking these dynamics into account is crucial for 
self-consistency of the calculation. 

In this paper the propagator for the graviton in quasi de Sitter space has been 
constructed in arbitrary dimension. In order to achieve this, first the propagator for a 
scalar field has been obtained. After a simple rescaling by a factor {ad)^^~^^'^^'^ , the 
scalar propagator in quasi de Sitter space has an identical structure as the 
corresponding scalar propagator in de Sitter space. The only difference is an order e 
correction to the effective mass. Due to this correction, the massless limit is regular, 
where it is singular in de Sitter space. 

Just as in de Sitter space, it is found that the graviton propagator in any FLRW 
space can be written as the sum of three scalar propagators with different amounts of 
coupling to the background geometry multiplying three different tensorial structures. 
If the background is quasi de Sitter, we are able to explicitly construct these scalar 
propagators. The graviton propagator in quasi de Sitter space allows for 
self-consistent calculations of graviton back-reaction and for quantum studies of 
tensor cosmological perturbations. 
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A Geometric quantities 



In this appendix we summarize some of the geometric quantities used in the text. 
The results in this appendix are valid in any FLRW spacetime. In this paper we work 
(mostly) with conformally flat space-times where the metric is given by 

g^Jiu = a^'n^y r]^y = diag(-l, 1, 1, 1) (108) 

We define the covariant derivative asfSH] 

V^F"/3 = d^y^f, + r-,y% - T'^^V\ (109) 

where the connection coefficients are given by 

r,^ = ^ ((^^5° + 5^5° + 5^7^,,) , (110) 

and a prime indicates a dervative with respect to conformal time. The Riemann 
tensor is defined as 

R'^^iiBu = d^v^^ - di^Tp^ + r^;^r^^, - r";^r^^. (iii) 

Using (jllOp we thus find 



' a' \ 2 / 



a 

RflV R flQU 



(112) 
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